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ABSTRACT 

We investigated the formation of arc-like structures in the infalling envelope around 
protostars, motivated by the recent Atacama Large Millimeter/Submillimeter Array 
(ALMA) observations of the high-density molecular cloud core, MC27/L1521F. We 
performed self-gravitational hydrodynamical numerical simulations with an adaptive 
mesh refinement code. A filamentary cloud with a 0.1 pc width fragments into cloud 
cores because of perturbations due to weak turbulence. The cloud core undergoes 
gravitational collapse to form multiple protostars, and gravitational torque from the 
orbiting protostars produces arc structures extending up to a 1000 AU scale. As well 
as on a spatial extent, the velocity ranges of the arc structures, ^ 0.5 km s -1 , are in 
agreement with the ALMA observations. We also found that circumstellar disks are 
often misaligned in triple system. The misalignment is caused by the tidal interac¬ 
tion between the protostars when they undergo close encounters because of a highly 
eccentric orbit of the tight binary pair. 

Key words: hydrodynamics - ISM: clouds - ISM: kinematics and dynamics - stars: 
formation - turbulence. 


1 INTRODUCTION 


It has been thought that multiple stars are formed through 
fragmentation from a high density portion of a molecular 
cloud core. The fragmentation occurs deep inside the molec¬ 
ular cloud c ore and it has not ever been observed direct ly. 

MC 27 (iM izuno et al. 19941 : 1 Onishi et al .111 9991 . 1 2002h or 
LI52IF (|Codella et al.ll 1997l l is a dense cloud core in Taurus 
and observations have suggested th at MC27/ L1521F is in a 
very e arly stage of s t ar for mation (|Onishi et a DE999). Re¬ 
cently iTokuda. et. al.l d2014h performed Atacama Large Mil- 
limeter/Submillimeter Array (ALMA) observations on dust 
continuum emission and molecular rotational lines toward 
MC27/L1521F. One of their new findings is an arc structure 
with a ~ 2000 AU length at the center of the molecular cloud 
core, suggesting that the arc structure can be attributed for 
a dynamical interaction between envelope gas and dense gas 
condensations. They also suggested that the turbulence pro¬ 
motes fragmentation of cloud cores during the collapse as 
show n by the numerica l simulations of turbul ent cloud cores 
(e.g., iBate et ahl 120021 : iGoodwin et aLlfeoQ^ j. In these nu¬ 
merical simulations, the orbital motion of fragments excite 
spiral arms in the center of the cloud core. However, these 
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studies focus mainly on the physical fragmentation process, 
while the observable features, such as the velocity structure 
of spiral arms, are still poorly understood. 

In this paper, the protostellar collapse of a turbulent 
core is calculated from a parent filamentary cloud by us¬ 
ing adaptive mesh refinement (AMR) simulations. We re¬ 
produced arc structures on a scale of about 1000 AU at the 
central part of the molecular cloud core. Based on these sim¬ 
ulations, the formation mechanism of the arc structures is 
discussed. This paper is organized as follows. In section [2j 
the models and simulation methods are presented. The re¬ 
sults of the simulations are shown in section [3j and they are 
discussed in section [H 


2 MODELS AND METHODS 

A filamentary molecular cloud is considered as the initial 
condition for this study because recent observations have 
revealed tha t filaments are basic compositions of molecular 
clouds fe.g.. lAndre et al1l20ld : lHacar et ai1l2013h . The fila¬ 
ment is assumed to be infinitely long and in an equilibrium 
state where thermal pressure supports an isothermal cloud 
against its self-gravity. The cloud therefore has a density dis¬ 
tribution of p(R) = po(l + R 2 /Ro)~ 2 , where R , po and Ro 
are a cylindrical radius, density on the filamentary axis, and 
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a scale height of th e filament, respectively (|Stodolkiewiczl 
1 1 963l : [Ost rikerll 1 964l i . We set R c = 0.05 p c, which mimics fila¬ 
ment s observed by the Herschel survey (| Arzoumanian et al.l 
1201 lh . The gas temperature was assumed to be T = 10 K 
at the initial stage and the corresponding sound speed is 
c s — 0.190 km s -1 . The density on the filamentary axis is 
given by po — 2 c 2 s /{'kGRq ) = 1.45 x 10 -19 gcm -3 (the cor¬ 
responding number density is no = 3.79 x 10 4 cm -3 ). The 
barotropic equation of state was assumed to be P(p) = 
c 2 p-\-Kp 7 / 5 with n = CgPci -2 ^ 5 , where the critical density is set 
at p cr = 10 -13 gem -3 (the corresponding number density is 
n cr fjt 2.6 2x 10 1 Q cm" 3 ), which was t a ken from th e numerical 
results of lMasunaga. Mivama. fc Inutsukal (| 19981 ). The mag¬ 
netic field was ignored for simplicity. The computational do¬ 
main is a cubic box with side lengths L = 4A ma x = 1.56 pc, 
where A max (= 7.81-Ro) is the wavelength of the most un- 
stable pert urbation against fragmentation of the filament 
(jNagasawal Il987h . Periodic boundary conditions were im¬ 
posed. 

Turbulence was imposed at the initial stage. The initial 
velocity field is incompressible wit h a power spectrum o f 
P(k) oc AT 4 , generated according to iDubinski et al] (Il995l b 
where k is the wavenumber. This power spectrum results 
in a velocity dispersion of < 7 (A) oc A^ 2 , in agreement with 
the Larson scaling relations (|Larsonlll98ll b The root mean 
square (rms) Mach number in the computational domain 
was set at unity. Note that the rms Mach number on a 0.1 pc 
scale is expected to be ~ 0.25 according to the scaling re¬ 
lations, and the turbulence is therefore subsonic on the fil¬ 
ament scale. Such a subsonic turbulence was suggested by 
the narrow molecu lar line widths of dense cores in Taurus 
(lOnishi et al.|[l998h . 

The evolution of the filamentary cloud was cal¬ 
culated using thr ee-dimensional AMR code, SFUMATO 
(|Matsumotol 120071 5 . The hydrodynamic scheme was modi¬ 
fied to have a third order of accuracy in space and a second 
order in time. The computational domain is resolved on a 
base grid of l = 0 with 256 3 cells. The maximum grid level 
was set at l = 11. The cell width is Ax m in = 0.613 AU on the 
finest grid of l = 11, compared with Ax max = 1.26 x 10 3 AU 
on the base grid of l = 0. The Jeans condition was employed 
as a refinement criterion; blocks are refined when the Jeans 
length is shorter than 8 times the cell width, i.e., A j < 8 Ax, 
where A j is the Jeans length fc.f.. iTruelove et ~ai1ll997l b 

The sink particle is introduced as a sub-grid model 
of protostars. The de tailed implementa t ions of the sink 
particles are shown in iMatsumoto et al.l ( 2015 1. The crit¬ 
ical density for sink particle formation is set at p S mk = 
1 x 10 -11 gem -3 (n s ink = 2.62 x 10 12 cm -3 ), and the sink 
radius is set at r s i n k = 4Ax m i n = 2.45 AU. 


3 RESULTS 

Figure [T| shows the column density distribution along the y- 
direction for four different spatial scales at t = 6.90 x 10 5 yr, 
or t p — 1.760 x 10 4 yr, where t p denotes an elapsed time af¬ 
ter the first sink particle formation. Figure [lja) shows the 
whole computational domain. The filamentary cloud frag¬ 
ments into the several cloud cores, one of which forms stars. 
The filamentary cloud is disturbed by the turbulence, which 
was imposed on the initial condition. Figure \Hb) shows the 


cloud core. The cloud core has a sharply peaked density dis¬ 
tribution in agreeme nt with molecular l ine observations to¬ 
ward MC27/L1521F (|Onishi et al.lll99~9h . The sink particles 
are associated with the peak position. The cloud core has a 
three-dimensional prolated shape an d it exhibits smooth iso- 
densitv surfaces, in agreement with Matsumoto & Hanawai 

m . 

On a smaller scale as shown in Figure [l][c), a triple sys¬ 
tem of sink particles can be seen. The sink particles are asso¬ 
ciated with several spiral arms, which correspond to the arc 
structures on a 1000 AU scale. The spiral arms are caused 
by the gravitational interaction between the sink particles 
and the envelope gas. The gas is accelerated by the gravi¬ 
tational torque from the orbiting sink particles, and the gas 
has a supersonic velocity of typically ~ 0.5 km s -1 , excit¬ 
ing shock waves of arc-shapes. Figure [2] shows the distri¬ 
bution of the line-of-sight velocity, which was evaluated as 
a first moment of velocity, f v y (r)p(r)dy/ J p(r)dy , assum¬ 
ing that the gas i s optically thin . This figure is comparable 
to Figure 4 from ITokuda et all ( 2014 b The arc shape and 
the velocity range ar e in agreement with the observations of 
ITokuda et alJ (2014). 

Note that the disk-like structure on a 1000 AU scale, 
which is shown as a region with log A [cm -2 ] > 23.5 in Fig¬ 
ure [U[c), is an infalling envelope. Its structure is not sup¬ 
ported by rotation, and the rotation velocit y is comparable 
to the infall velocity there. The Q-value of iToomref (Il964h 
is larger than unity, indicating that the infalling envelope as 
well as the arcs are not significantly self-gravitating. 

Figure [ljd) is a close-up view of the triple system. The 
sink particles are labeled by identification numbers 0, 1 and 2 
in order of formation epoch (hereafter referred to as Sinks 0, 

1, and 2). Sinks 0 and 2 constitute a close binary with a 
separation of 20 — 30 AU, while Sink 1 orbits around the 
close binary with a long separation of 100 — 200 AU. At the 
stage as shown in Figure [lj Sinks 0, 1, and 2 have masses of 
0.20 M©, 0.32 M©, and 0.13 M©, respectively. Sink 1 has the 
largest mass among the sink particles, and it has the largest 
circumstellar disk with a radius of ~ 30 AU, while the other 
sink particles have small circumstellar disks with radii of ~ 
7 AU. The circumstellar disk of Sink 1 is relatively massive; 
the mass ratio between the disk and the sink particle remains 
at Mdisk/Mstar ~ 0.2 for a long period of t p > 700 yr, while 
other sink particles take ratios of Mdi s k/M s t a r ~ 0.003 — 0.01. 
The massive disk of Sink 1 has the spiral arms caused by the 
gravitational instability; the disk exhibits the Q-value less 
than unity. They promote a high accretion rate onto Sink 1, 
- 2 x 10 -5 M©yr -1 . 

The circumstellar disk of Sink 2 is inclined with respect 
to the other circumstellar disks. The misalignment of the 
circumstellar disk is caused by tidal interactions from Sink 0 
whe n Sink 2 unde rgoes sccessive close encounters with Sink 0 
fc.f.. iHellerl1 1993l b These close encounters are caused by the 
highly eccentric orbit of the tight binary pair of Sinks 0 and 

2. During the evolution, rapid changes in the disk orientation 
occur several times for the disks of Sinks 0 and 2. 

Figure [3] shows the changes in separation between all 
of the pairs of sink particles. In the early stages, four sink 
particles were formed and they exhibit chaotic orbits (Fig¬ 
ure [3][b)). In their chaotic orbital motion, Sinks 1 and 3 
merge at t p = 3800 yr (the cyan line). Subsequently, Sinks 1 
and 2 form a close binary with a separation of 10 — 40 AU 
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Figure 1. Column density distribution along the y-direction on four different spatial scales at t = 0.690 Myr ( t p = 1.760 X 10 4 yr). 
The color scales depict the column densities. The black crosses indicate the positions of the sink particles. In panels (c) and (d), the 
coordinates are offset so that the center of each panel coincides with the center of mass of the sink particles. The sink particles are labeled 
with identification numbers 0, 1, and 2 in panel (d). 


(the blue line). At t p = 6000 yr, Sink 0 approached the close 
binary of Sinks 1 and 2. Sink 0 and 1 were then exchanged 
to from a new close binary with Sinks 0 and 2. The orbit of 
Sink 1 then gradually decreases the eccentricity from 0.3 to 
0.2 (the blue and red lines in Figure [3]( a)). The close binary 
of Sinks 0 and 2 has a high eccentricity (the green line), 
leading to the misaligned circumstellar disks. 

Because the sink particles have a chaotic epoch after 
their formation, the density structure in the infalling en¬ 
velope, such as the arc, depends on the realization of the 
initial turbulence. We simulated the same model as the pre¬ 
vious one but changed the random seed when generating 
the turbulence. Figure H] shows the column density distri¬ 
bution at t = 8.15 x 10 5 yr (t p = 2.14 x 10 4 yr). By this 


stage, this realization model has produced six sink particles 
in total and they eventually merge into three sink parti¬ 
cles. The sink particles undergo chaotic orbital motions until 
t p ~ 2 x 10 4 yr, which is a considerably longer period than in 
the previous model. This realization model also exhibits arc 
structures on a 1000 AU scale. The line-of-sight velocities 
of the arc structures also range within ~ ±0.5 km s -1 . This 
model also exhibits misaligned circumstellar disks caused by 
close encounters between sink particles. In the stage shown 
in Figure [H the most distant sink particle (Sink 0) has a 
circumstellar disk inclined with respect to the other circum¬ 
stellar disks. 
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Figure 2. Line-of-sight velocity distribution in the arc struc¬ 
tures. The color scale shows the velocity along the line of sight. 
The contour shows the column density on a logarithmic scale with 
log(A/'/cm -2 ) = 23.4, 23.6, • • •, where N denotes the column den¬ 
sity. The velocity is shown only in the regions where the column 
density is greater than the lowest contour level. 



4 SUMMARY AND DISCUSSION 

We demonstrated the formation of arc structure on a 
1000 AU scale in the infalling envelope during the early 
phase of star formation as a consequence of the collapse 
of a weakly turbulent cloud core. The density and velocity 
structures reproduced here are in agreement with the arc 
structure revealed by the recent ALMA observations toward 
MC27/L1521F. We also found that the misaligned disks are 
formed in triple s ystems. _ 

According to iTokuda et al.l (|2014h . MC27/L1521F has 
three high-density condensations observed by the dust con¬ 
tinuum, one of which ( MMS-1) was also observed by the 
Spitzer space telescope (jBourke et al.ll2006h . The other two 
condensations (MMS-2 and MMS-3) have high densities of 
10 6-7 cm -3 , and MMS-2 is a candidate for the first core as 
predicted bv lLarso 3 (Il969l) . 

A possible interpretation for these high-density conden¬ 
sations is that they are formed by fragmentation during the 
protostellar collapse of the cloud core as shown in our simu¬ 
lations. The velocity of the arc structure is ranged typically 
within ~ ±0.5kms _1 , and it is consistent with the dynami¬ 
cal velocity of the system, e.g., the Kepler velocity for a mass 
of 0.2 M© and a radius of 500 AU. Moreover, the time scale of 
the arc structure is 1000 AU/0.5 kms -1 ~ 10 4 yr, in agree¬ 
ment with the epoch shown in Figure [I] t p = 1.760 x 10 4 yr. 
However, a few issues remain unsolved in the following 
points. First, the masses of the sink particles grow exceed¬ 
ing 0.1 Mq with this t ime scale . The se ma sses are greater 
than those reported bv ITokuda et al.l (|2014h . Second, MMS- 
2, the first core candidate, may have a short lifetime of less 
than 1000 yr a ccording to the recen t radiation-transfer MHD 
simulations bv lTomida et al.l (l2013h . Within such a short pe¬ 




Figure 3. Separations between all of the pairs of sink particles as 
a function of time after the first sink particle formation. The label 
n — m associated with each line indicates the separation between 
Sinks n and m. Panel ( b ) is an enlargement of the chaotic period 
in the early phase. 


riod, an arc structure can not extend to a 1000 AU scale with 
its typical velocity of ~ 0.5 km s -1 . 

Protostars undergo chaotic orbits in the period of sev¬ 
eral xlO 3 — 10 4 yr after their formation. Their chaotic or¬ 
bits influence the evolution of multiple stars as well as the 
surrounding envelope, such as the arc structures. The grav¬ 
itational few-body problem is therefore a key issue in the 
early phase of low-mass star formation. This is somewhat 
analogous to the interactions in “the competitive accretio n 
scenario” for massive star formation (|Bonnell et al.ll2004i L 
in which the gravitational N-body problem is supposed to 
govern the dynamics of clusters and therefore accretion. 

Another scenario to consider for the formation of an 
arc structure is a magnetic wall around the protostars. The 


formation of a ma 
retical works (e.g., 

gnetic wall was predicted by some theo- 


z McKed 19961: iTassis & Mouschovias 

200fil: Izhao et al. 2011 

). The three-dimensional simulations 


showed that the magnetic wall produces cavities in the in¬ 
falling envelope. The rims of the cavities have higher den¬ 
sities than the envelope and they may be observed as arc 
structures. However, the magnetic field suppresses fragmen- 
tation for the typical parameters of molecular cloud cores 
(Ma chida et al.l l2005h . The formation mechanism of the 
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Figure 4. Column density distribution showing arc structures. 
The black crosses indicate the positions of the sink particles, 
which are labeled with identification numbers. The realization 
is different from that shown in Figure [l] 


high-density condensations, MMS-2 and MMS-3, can then 
not be explained by the fragmentation process. The observed 
arc structure therefore implies that the magnetic field is 
weak enough to promote fragmentation of the cloud core. 

Misalignment of a circumstellar disk is caused by the 
tidal interaction with sink particles when they undergo 
close encounters. The close encounters are promoted by the 
chaotic orbits in a triple system. Thus, triple systems are 
expected to have misaligned circumstellar disks, and to also 
be associated with arc structures if they are young enough 
to be embedded in the dense envelope. 

Misaligne d disks were reported for the youn g binary sys¬ 
tem, HK Tau (I Jensen V Akesonll2014lh AS 205 (jSalvk et al J 
l2014h . and V2434Qri (I Williams et al .1120141 5 based on the re¬ 
cent ALMA observations. The young triple system, T Tau, 
includes the tight binary of T Tau Sa and T Tau Sb, and the 
existence of a misaligned circumstellar disk has been indi¬ 
cated therein (ISkemer et al.l 120011 iRatzka et al.ll2009h . Our 
model suggests that this misaligned disk could be caused 
by the tidal interaction between the tight binary pair in the 
triple system. 
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